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Abstract 

Without an ultraviolet cut-off, the time evolution of the classical Yang-Mills equations 
give rise to a never ending cascading of the modes towards the ultraviolet, and ergodic 
measures and dynamical averages, such as the spectrum of characteristic Lyapunov 
exponents (measures of temporal chaos) or spatial correlation functions, are ill de- 
fined. A lattice regularization (in space) provides an ultraviolet cut-ofi' of the classical 
Yang-Mills theory, giving a possibility for the existence of ergodic measures and dy- 
namical averages. We analyze in this investigation in particular the scaling behavior 
(5 = d log A/d log S of the principal Lyapunov exponent with the energy of the lattice 
system. A large body of recent literature claims a linear scaling relationship (/? = 1) 
between the principal Lyapunov exponent and the average energy per lattice plaquette 
for the continuum limit of the lattice Yang- Mills equations. We question this result 
by providing rigorous upper bounds on the Lyapunov exponent for all energies, hence 
giving a non-positive exponent, /? < 0, asymptotically for high energies, and we give 
plausible arguments for a scaling exponent close to /3 ~ 1/4 for low energies. We argue 
that the region of low energy is the region which comes closest to what could be termed 
a "continuum limit" for the classical lattice system. 
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The classical Yang-Mills equations for non-Abelian gauge fields provide an interesting 
class of dynamical systems for which non-linear self-coupling terms open up the possibility 
of chaotic behavior even in toy-models where the gauge fields are spatially homogeneous 
[01 " 0- However, a qualitative important new dynamical feature comes into play when 
one considers the spatially inhomogeneous classical Yang-Mills equations: There is a never 
ending cascading of the dynamical degrees of freedom towards the ultraviolet, generated by 
the time evolution of the Yang-Mills equations as we shall discuss below. In spite of this 
"ultraviolet catastrophe" , we note that the solutions are well behaved, in the sense that there 
are no "finite time blow up of singularities". For example, the class of classical Yang- Mills 
equations in 3+1 dimensions is known to have solutions in suitable Sobolev spaces, such that 
their norms do not blow up to infinity in finite time 0, H]. Without an ultraviolet cut-off of 
the Yang-Mills field equations, it is however not possible - even when the fields are confined 
to a finite volume - to define equilibrium statistical mechanics, i.e. make reasonable sense out 
of infinite-time time averages of observables and correlation functions. The reason for this is 
that although the system has finite space volume, the phase space volume is still infinite. As 
is well-known in the case of electromagnetism, a reasonable definition of equilibrium leads 
to an ultraviolet catastrophe as seen in the Rayleigh- Jeans law for energy (per frequency) 
density. The arguments in the electromagnetic case are purely thermodynamic and do not 
involve the classical evolution of the fields. In fact, since the corresponding equations are 
linear the problems would not show up dynamically unless one couples the fields to charged 
particles. In contrast, the Yang-Mills equations for non-Abelian gauge fields include non- 
linear self-coupling terms which open up the possibility for a chaotic behavior in the classical 
evolution and in the non-homogeneous case it leads to the infinite cascade of energy from the 
long wavelength modes towards the ultraviolet. This is a natural dynamical interpretation 
of the ultraviolet catastrophe if one assumes that on average a trajectory should simulate a 
micro canonical ensemble where, in an equilibrium situation, all modes have the same energy 
on average. But due to the infinite phase space volume, this implies that this average energy 
per mode has to be zero (!), whence the cascade. We would like to add that this picture is 
well supported by the articles of Furusawa |0 and Wellner [|l^ . 

This tendency of the mode frequencies cascading towards the ultraviolet will completely 
dominate the qualitative behavior of the classical Yang-Mills equations, and the "ultraviolet 
catastrophe" has for some time been emphasized by us (H.B.N, and S.E.R., cf. e.g. dis- 
cussion in |26|) as a major obstacle to simulate the classical continuum Yang-Mills fields in 
a numerical experiment over a long time span. (This obstacle has in our opinion received 
insufficient attention in the various studies attempting at discussing and modeling chaotic 
properties of spatially inhomogeneous classical Yang-Mills fields). There is no mechanism, 
within the classical equations, which prevents this never ending cascading of the modes 
towards the ultraviolet. Nature needs h, the Planck constant, as an ultraviolet regulator. 
Indeed, both Abelian and non-Abelian gauge fields are implemented as quantum theories in 
Nature. 

Here, however, we shall discuss the possibility of using a lattice cutoff in a purely classical 
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treatment to regularize the equations and giving sensible results in the limit when the lattice 
spacing a tends to zero. The model we consider is the Hamiltonian lattice formulation first 
introduced by Kogut and Susskind (in a quantum context). There is still a cascading of 
modes towards the ultraviolet, i.e. towards the lattice cut-off, and this ultraviolet cascade 
will still dominate the dynamical evolution of smooth initial field configurations. However, 
in this Hamiltonian formulation on a large but finite lattice, the phase space is compact for 
any given energy and thus the system can reach an equilibrium state among the modes (a 
'thermodynamic equilibrium'). The lattice regularization of the theory opens up for defining 
dynamic and thermodynamic properties, which are not defined in the classical Yang-Mills 
field theory without regularization. It could e.g. be ergodic (modulo constraints) with re- 
spect to the Liouville measure, in which case it makes sense to talk about its micro-canonical 
distribution and approximating this by looking at 'typical' classical trajectories. The fun- 
damental assumption of thermodynamics asserts that on average the two approaches give 
the same result if we have a large system, and we may then naturally introduce correlation 
functions and possibly a correlation length ^ (measured in lattice units) of the system. One 
hopes to define a continuum theory if, by judicious choice of the parameters in the system, 
one obtains a physical correlation length in the limit when the lattice constant goes to zero, 
i.e. 

^(a,^(a),...) X a ^ £ ^ , as a ^ 0. (1) 

In equation (|1|) the correlation length ^ in the lattice system is a function of lattice model 
parameters such as lattice spacing a, average energy density -E(a), etc. Condition (|I]) implies 
that the correlation length diverges when measured in lattice units and only if this is the 
case do we expect the lattice system to lose its memory of the underlying lattice structure. 
Condition (|I|) is, of course, just a necessary, not sufficient condition for the lattice system 
to approach the continuum theory. For example, the limiting system might still carry the 
symmetries of the lattice. We shall return in more detail to the important issue of what we 
could mean by a "continuum limit" of a classical lattice regularized non-Abelian gauge field 
theory. 

The Kogut- Susskind Hamiltonian is obtained from the Wilson action, approximating the 
Yang-Mills equations, by keeping a lattice grid in space while letting time become continuous. 
For a derivation we refer to |^ and also to [|T^ from which we adapt our notation. We will 
restrict attention to the lattice gauge theory in 3-1-1 dimensions based on the gauge group 
SU{2). One considers a finite size 3 dimensional discretized box, i.e. a finite lattice having 
A^^ points where nearest neighbor points are separated by a distance a > 0. The phase space 
corresponding to this is then a fibered space where the tangent manifold of the Lie-group 
SU{2) is assigned to each of the links, i E A, connecting nearest neighbor lattice points. 
More precisely, we have to each i G A associated a link variable Ui G SU{2) as well as its 
canonical momentum Pj G Tjj^SU{2). A point in the entire phase space will be denoted 

X = {f/„ P,},eA G M = \{T SU{2) . (2) 

A 
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The resulting Kogut-Susskind Hamiltonian generating the time evolution of the orbit X(i:) 
can be written in the following way, where for simplicity we omit the coupling constant factor 
2/ g"^ which anyway is arbitrary in a classical theory : 

H{a, X('^)) = - 5: ^tr(P.i^t) + i - LtUu) • (3) 

Here the last sum is over elementary plaquettes bounded by 4 links and Un denotes the 
path-ordered product of the 4 gauge elements along the boundary of the plaquette □. The 
last term, the potential term, is automatically bounded and for a given finite total energy 
the same is the case for the first term, the kinetic term. Thus the phase space corresponding 
to a given energy-surface is compact. 

The compactness of the phase space implies that the spectrum of Lyapunov exponents 
(which we overall will assume to be well defined quantities for the lattice system) is indepen- 
dent of the choice of norm on the space of field configurations. The measures of distances 
between lattice field configurations X(t) and X(t) which are employed in - |19| are 
semi-positive and, typically, of the form (cf. e.g. fll 



VO^X) {Y.MP,Ph -tr{P:P})\ + Y.\t^Un-tTUu\) (4) 

which in the limit as the lattice constant a — > measures the (average) local differences in 
the electric and magnetic field energy, 

V{{E,B),{E,B))^^ j d?x{\E'-E'\ + \B'-B'\) . (5) 

Gauge equivalent field configurations have a vanishing distance since the distance measure 
(H) is gauge invariant .Q As discussed at the end of the appendix such a choice of distance 
will not affect the calculation of the largest Lyapunov exponent. 

It is easy to see that if X*^"^^) (t) solves the Hamiltonian equations for a = 1 then X*^"^) (t) = 
X('^=i)(t/a) solves the same equations for general a and this time-scaling implies the following 
scaling of the maximal Lyapunov exponent :0 

1 J|5X(°) 
Amax(a) = hm - log 



^^t '''115X^(0)11' 
1 1 J|5X("=^(t/a)|| ^ X^a^ja = 1) 
a W^) ^"g^ PX(--^)(0)|| ^ = a 



^As regards gauge invariance in the resulting measures of chaos, we note that since "space-time" is not 
involved in the gauge transformations, we expect Lyapunov exponents and correlation lengths to be gauge 
invariant measures of temporal and spatial chaos in the context of Yang-Mills fields, whereas such gauge 
invariant measures of chaos are more difficult to construct for a typical metric in general relativity ||27| . 

^In the equation (^) we should in fact replace the small deviation vector (5X by a tangent vector to the 
manifold, cf. the appendix. 



4 



Equation implies that aX^axi^i) is invariant under rescaling of a. If we assume ergodicity 
over the energy surface (modulo the known constraints) we see that in particular, a\max{.a) 
is uniquely determined by the energy aH{a) = H{a = 1) or equivalently when lattice size 
A^^ is fixed it has to be some function of the energy per plaquette aE{a) = E{a = 1), 

a Xmax{a, E{a)) = f {a E{a)) (7) 

Such a functional relationship even persists in the thermodynamic limit, N oo, for any 
fixed values of a and E{a). Thus in order to study the dependence of the maximal Lya- 
punov exponent with energy density, it is sufficient to consider the equations of motion 
for a fixed value of the lattice constant a, e.g. a = 1, as a function of energy density (oc 
energy/plaquette) and rescale the results back afterwards. 

Studying the classical dynamics in time rather than the ensemble distribution makes it 
possible to extract dynamical information as well as thermodynamic information. By mea- 
suring the maximal Lyapunov exponent Amax(a) for some fixed lattice spacing a of the system 
as a function of the energy, one may also hope to gain insight into the spatial correlations of 
the fields. Although a precise relationship has not been established yet, it seems plausible 
that when propagation of information to nearest neighbors occurs at a fixed speed, a small 
value of Xmax should lead to a larger correlation length than a larger value of Xmax- The rea- 
son is that in order for the dynamics to create long range spatial correlations in the system, 
information has to propagate for a long time (= distance) without seriously being attenuated 
by the chaotic behavior, the 'strength' of which is refiected in the value of Xmax- Clearly, it 
deserves further investigation to establish a precise relationship between Xmax and the spatial 
correlation length ^ (in lattice units). Such investigations would be of considerable interest 
also in the more general case where the non-Abelian gauge fields are given a non-zero mass, 
with an adjustable mass parameter, by coupling to a Higgs-field. In our case, which involves 
massless non-Abelian gauge fields, we expect that it is in the limit as Xmax goes to zero that 
the spatial correlation length ^ (in lattice units) will diverge. f\ 

In fact, a complicated question (which needs further study) is what we could mean by a 
"continuum limit" for a classical, lattice regularized non-Abelian gauge field theory. Whereas 
the issue of extracting expectation values of observables in the "continuum limit" is stan- 
dard in interpreting simulations of quantum field theories in the imaginary time (Euclidean) 
formalism, the issue of extracting "continuum results" (rather than lattice artifacts) is much 
less well-known in the case of real-time simulations of classical non-Abelian gauge theories. 
Moreover, it is an important question of interest in physics, since not many non-perturbative 

^It is a generic expectation for a chaotic, non-integrable, spatially extended system (with a spatial propa- 
gation of disturbances at some fixed speed) that increasing the temporal chaos, as captured e.g. by measuring 
the maximal Lyapunov exponent Xmax, will be connected to more spatial chaos, as measured, say, by an 
inverse spatial correlation length ~ Cf. e.g. discussions in | 
spatially extended Hamiltonian system can also have finite spatial correlations due to the initial conditions 
rather than the dynamics. I.e., an integrable Hamiltonian spatially extended system will not be able to 
establish spatial correlations in field configurations which initially are spatially de-correlated. 



24, pa. Note, however, that an integrable, 
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methods are available to study the time evolution of quantum field theories, such as Q.C.D., 
or the weak sector of the Standard Model, and one almost has to resolve to study the real 
time evolution of the classical fields - an approach approximately justified when probing 
physical situations (e.g. at high temperatures) where the quantum fields are expected to 
behave semi-classically (see also, e.g. discussions in 

For the study of the time evolution of Yang-Mills fields which initially are far from 
an equilibrium situation, the (classical) field modes will exhibit a never ending dynamical 
cascade towards the ultraviolet and after a certain transient time, the cut-off provided by 
the spatial lattice will prevent the lattice gauge theory from simulating this cascade. It is 
therefore immediately clear that the lattice regularized, classical fields will not approach a 
"continuum limit" in the sense of simulating the dynamical behavior of the classical contin- 
uum fields in the t ^ oo limit (we may also say that the continuum classical theory does 
not exist in the infinite time limit). 

We have several different forms of lattice artifacts in the lattice simulation of real-time 
dynamical behavior of the continuum classical Yang-Mills fields: 

(1) Lattice artifacts due to the compactness of the group. The magnetic term (the second 
term) in the Kogut-Susskind Hamiltonian (|^) is uniformly bounded, < 1 — ^Trlln < 2, due 
to the SU (2) compactification. Thinking in terms of statistical mechanics for our classical 
lattice system, we expect that after some time the typical field configuration has equally much 
energy in all modes of vibration - independent of the frequency^ - and the total amplitude of 
the classical field, and the energy per lattice plaquette, is thus small for a fixed low energy. 
For low energy, when the average energy per plaquette is small, the lattice artifacts due to 
the compactness of the gauge group are thus negligible. 

(2) For small energy per plaquette, we thus expect that the dominant form for lattice 
artifact is due to the fact that an appreciable amount of the activity (for example the energy) 
is in the field modes with wavelengths comparable to the lattice constant a. This short 
wavelength activity at lattice cut-off scales is unavoidable in the limit of long time simulation 
of an initially smooth field configuration (relative to the lattice spacing), or already after a 
short time if we initially have an irregular field configuration. 

We conclude that for the simulation of gauge fields far from equilibrium, we will have the 
best "continuum limit" if we simulate, for a short period of time, an initial smooth ansatz 
for the fields in the region of low energy per plaquette. Q 

^This situation is very different in the quantum case. Planck's constant h introduces a relation E = hco 
between the energy of a mode of vibration and its frequency, implying that a mode with a high frequency 
also has a high energy. With a given available finite total energy, modes with high frequencies will therefore 
be suppressed. Quantum mechanically, we thus have that at low energy only excitations of the longest 
wavelengths appear. 

^We are here imagining a situation where the spatial correlations in the monitored field variables are 
so large that they lose memory of the underlying lattice structure (including the lattice spacing a). In the 
extreme opposite limit, one could imagine models in situations with random fluctuating fields on the scale 
of the lattice constant, i.e. with (almost) no spatial correlations from link to link. If the field variables 
fluctuate independent of each other (independent of their neighbors) one could imagine the model to be 
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This is if we have no quantum mechanics. However, Yang-Mills fields exist as quantum 
theories. Thus, the interesting physical definition of a "continuum limit" of the classical 
lattice gauge theory, is to identify regions in the parameter space for the classical lattice 
gauge theory which, for (shorter or longer) intervals of time, probe the behavior of the time 
evolution of semi-classical initial configurations (with many quanta) of the quantum theory, 
for example Q.C.D. (See also e.g. discussion in [p^ ). 

There is by now a large amount of literature which reports investigations of temporal 
chaos on the lattice gauge theory measured by the spectrum of Lyapunov exponents (either 
the principal Lyapunov exponent or the entire spectrum of Lyapunov exponents) and how 
Lyapunov exponents depend on the average energy per plaquette, E, of the lattice field 
theory. We refer to a sequence of articles by Miiller et al fl^. Gong [l^. Biro et al. |lj 



and a recent book by Biro et al [r^| which present numerical analyses of the classical SU{2) 



lattice gauge model. Their numerical results provide evidence that the maximal Lyapunov 
exponent is a monotonically increasing continuous function of the scale free energy/plaquette 
with the value zero at zero energy. 

Miiller et al. report a particular interesting interpretation of numerical results for 
the dynamics on the lattice, namely that there is a linear scaling relation between the scale 
free maximal Lyapunov exponent, Xmax{ci = 1) and the energy per plaquette E{a = 1). The 
possible physical relevance of this result is seen when we rescale back to a variable lattice 
spacing a and note that the observed relationship is in fact a graph of aXmax (o) as a function 
of aE{a). Thus being linear, cancellation of a factor a implies that Xmax{0') = const x E{a) 
and thus there is a continuum limit a — >■ either of both sides simultaneously or of none of 
them. In the particular case where the energy per mode (oc energy per plaquette) is taken 
to be a fixed temperature T (cf. [^), one deduces that the maximal Lyapunov exponent 
has a continuum limit in real time. 

We shall here, however, note that the linear relationship found between Xmax and E is 
based on a graph of their interdependence (indicating that such a scaling relation exists) in 
the approximate scaling region (for a = 1/2) of E between 0.6 and 4, where the maximal 
Lyapunov exponent grows from 0.2 to 1.3 (cf. figure 1) : 

X^ax ~ 0.32 E . (8) 

We shall argue that the apparent linear scaling relation is a transient phenomena residing 
in a region extending at most a decade between two scaling regions, namely for small energies 
where the Lyapunov exponent scales with an exponent which could be close to 1/4 and a 
high energy region where the scaling exponent is at most zero: 

^ dlogX ^ r -1/4 for E ^ , . 

^ dlogE 1 <0 for^^oo ^ ^ 



invariant (with respect to the monitoring of many variables) under changes of the lattice spacing, a. Thus, it 
appears that lattice cut-off independence of numerical results can not be a sufficient criterion for the results 
to report "continuum physics" . 
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In fact, we believe that the apparent hnear scahng region, equation (^, is a complete arti- 
fact of the SU{2) compactification on a lattice. We argue by giving a newly obtained ||2^ 



rigorous upper bound for the maximal Lyapunov exponent, as well as plausible lower bounds 
obtained from general scaling arguments ^6] of the continuum classical Yang-Mills equations 



in accordance with simulations on homogeneous models M and consistent with the figures in 



1^ and JT^. Thus the existence of the two scaling regions, equation (^, is established by a 



combination of numerical evidence and analytical arguments. In figure 1 (a) we have plotted 



the numerical results reported in Miiller et al [0] which exhibits a systematic deviation 
from a straight line. A recent numerical investigation by Krasnitz (cf. [^, figure 3) has 
independently confirmed such a systematic deviation. According to Krasnitz statistically 
significant deviations from the straight line is as large as 10 %. In figure 1 (b) we have made 



a log-log plot of the results obtained by Miiller et al |T| 



Regarding the limit for high energies per plaquette a rigorous result (a sketch is 
shown in Appendix A) shows that the SU (2) scale free (a = 1) lattice Hamiltonian in d 
spatial dimensions has an upper bound for the maximal Lyapunov exponent 



Ab = V(^^-l)(4 + V17) (10) 

which for d = 3 becomes Xb = 4.03... This result is arrived at by constructing an appropriate 
norm on the phase space and showing that the time derivative of this norm can be bounded 
by a constant times the norm itself, hence giving us an upper bound as to how exponentially 
fast the norm can grow in time. The upper bound ( [T0| ) shows that a linear scaling region, 
i.e. a constant j3 = dlogX/dlogE, can not extend further than around ~ 10 on the 
figure 1. Beyond that point the maximal Lyapunov exponent either saturates and scales 
with energy with an exponent which approaches zero or it may even decrease over a region 
of high energies, yielding a negative exponent^. It should be noted that the upper bound 
( p!OD is independent of the lattice size and the energy (but scales with I /a). 

In the opposite limit, the "continuum limit", where the average energy per plaquette 
E ^ 0, the finite size of the lattice makes it much more difficult to analyze the behavior 
of the gauge fields and the principal Lyapunov exponent on the lattice. However, we shall 
argue that the scaling exponent of the Lyapunov exponent more likely will be closer to ~ 1/4 
than the scaling exponent ~ 1 observed in the intermediate energy region 1/2 <E <4. 

In figure 1 (b) we have made a log- log plot of the results obtained by Miiller et al. |T3]. It 



appears that for ~ 1/2 there is a cross-over to another scaling region. Although the data 
points in this region are determined with some numerical uncertainty we note that they 
are consistent with a scaling with exponent ~ 1/4. Moreover, it is a numerically established 
fact that the homogeneous Yang-Mills equations have a non-zero Lyapunov exponent which 
by elementary scaling arguments scales with the fourth root of the energy density and has 
the approximate form [Q: 

A„,ax ~ 0.38 E^/^ (11) 
^ Some very heuristic arguments yields f3 = —1/2 asymptotically. 
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On the lattice the above scahng relation is valid for spatially homogeneous fields, i.e. the 
maximal Lyapunov exponent scales with the fourth root of the energy per plaquette. By 
continuity, fields which are almost homogeneous]] on the lattice will, in their transient, initial 
dynamical behavior, exhibit a scaling exponent close to 1/4. In fact, the same scaling 
exponent^ would also hold for the inhomogeneous Yang-Mills equations had the fields 
been smooth relative to the lattice scale, so derivatives, d^, are well approximated by their 
lattice equivalent and we are allowed to scale lengths as well. This is seen from scaling 
arguments for the continuum Yang-Mills equations : 

Not taking boundary conditions into account, these equations are invariant when and 
are scaled with the same factor a. That is, if A{x,t) is a solution to the equations, then 
^A{ax, at) is also a solution. The energy density E, which is quadratic in the Yang-Mills 
field curvature tensor F'^'^, then scales with a^. The same arguments as leading to equation 
then show that if we perform a measurement of the maximal Lyapunov exponent Xmax 
over a time short enough for the solutions to stay smooth, then Xmax scales with E^^^, as was 
also observed by Miiller et al. p. 3389 in |T^. These scaling arguments do not carry over to 



infinite time averages since solutions on the lattice tend to be irregular. We believe, however, 
that the cross-over observed in figure 1(b) is a feature of the Kogut-Susskind model and not 
just a finite size effect, in particular since there is evidence that the correlation length for 



energies ~ 1/2 is of the order of a few lattice units only. As the numerical simulations of 
Miiller et al. are performed on lattices ~ 20'^ and since the calculation of Lyapunov exponents 
is a local calculation (when considering infinitesimal deviations), a finite size effect should 
show up for energies somewhat below E ~ 1/2. 

The expectations, arrived at here, are supported by an even more striking result, which 
is the numerical simulations of Biro et al (fl^, figure 12) for the Kogut-Susskind lattice 
model with a U{1) group showing (see figure 2) a steep increase of the maximal Lyapunov 
exponent with energy/plaquette in the interval 1 < -E" < 4. The continuum theory here 
corresponds to the classical electromagnetic fields which have no self-interaction and thus 
the Lyapunov exponent in this limit should vanish. The discrepancies in this case were in 
1^ attributed to a combined effect of the discreteness of the lattice and the compactness of 



the gauge group U{1) and were not connected with finite size effects. This suggests strongly 
that in the case of SU{2) a similar cross-over around E 1/2 should have the same origin 
and likewise, not be seen as a finite size effect. In particular, it makes it difficult to believe 
that we in the case of numerical studies of a SU{2) Kogut-Susskind Hamiltonian system can 
base continuum physics on results from simulations in the same interval of energies where the 

''Note that such field configurations are (ungeneric) examples of field configurations which exhibit corre- 
lation lengths much larger than the lattice spacing. 

*In the lattice simulations the energy of a given plaquette is a fluctuating quantity and since the function 
E^^^ is convex one would expect the coefficient of proportionality to be somewhat smaller than in formula 
(O). Numerically the coefficient turns out to be around half of the above value, cf. figure 1 (a-b). 
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f/(l) simulations fail to display continuum physics. On the contrary, we suspect that what 
could reasonably be called "continuum physics" has to be extracted from investigations of 
the Kogut-Susskind lattice simulations for energies per plaquette which are at least smaller 
than E 1/2. In conclusion, we believe that there is indeed a cross-over around ~ 1/2 but 
we do not attribute this to finite size effects^ which would invalidate the physical relevance 
of the lower region. Rather, we believe and conjecture that this phenomena will persist 
even for infinitely big lattices still showing a cross-over around E ~ 1/2, below which we 
approach "continuum physics" and above which discreteness of the lattice and compactness 
of the gauge group are no longer negligible. 

There is quite a simple intuitive explanation for the kind of transition taking place for 
energies in the region above 1/2 for the SU{2) lattice gauge model and for the saturation 
of the Lyapunov exponent in the regime for high energy per plaquette: The Hamiltonian 
@ consists of two terms, the first one is a kinetic energy term (electric term) which is a 
quadratic form in the momenta, the second is a potential term (magnetic term) which is 
uniformly bounded < 1 — ^TrUn < 2, due to the SU{2) compactification. For small 
energies per plaquette, <^ 1, both terms will be of comparable size. In fact, as energy 
increases from zero, the potential energy increases too, but as long as the energy/plaquette 
is much smaller than 1, the gauge holonomy Uu calculated around a Wilson loop is close to 
the identity so that field curvatures are small. The dynamics is then essentially confined to 
the Lie algebra, i.e. the tangent plane of the identity element in SU{2). However, since the 
Lie-group itself is curved, higher order non-linear terms become important as the energy per 
plaquette increases further. Although this does not imply 'more chaos', it does suggest a 
steeper increase in the Lyapunov exponent. As energy increases even more, E ^ 1, the Un 
will experience the finiteness of the group SU{2) and energy will only be pumped into the 
kinetic term. We note that if only the kinetic term had been present in the Kogut-Susskind 
Hamiltonian (|^), the system would be integrable. Since the potential term only provides 
a uniformly bounded perturbation of the dynamics, it is reasonable to expect (and shown 
in the appendix) that the spectrum of Lyapunov exponents will saturate as the energy per 
plaquette increases for the Kogut-Susskind model. The chaos generated by the non-linear 
potential energy term does not increase, only the energy in the kinetic energy (the electric 
fields) increases. 

Finally, we will give a few comments on possible implications of the scaling relation (^ 
proposed in this investigation. Since by intrinsic scaling arguments (cf. equation (|^)) one 
has a functional relation between aX^axi^i) and aE{a), a 1/4 scaling for small energies would 
imply that A (a) oc a' -^/^E{aY/\ in which case one cannot achieve a continuum limit simul- 
taneously for the maximal Lyapunov exponent and the temperature (assuming that it is 
proportional to the average energy per plaquette E{a), cf. |^), in particular, the former 



would be divergent if the temperature is kept fixed. There is no particular contradiction in 



^It would be valuable to study the spatial correlation lengths on the lattice system for energies per 
plaquette in the interval < E< 1, i.e. for an interval of energies which includes the apparent cross-over at 
E - 1/2. 
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this statement, however, as there is, a priori, no need for having a finite Lyapunov exponent 
in the continuum hmit. The erratic and fiuctuating behavior of the fields one expects in 
time as well as in space (for numerical evidence, cf. also [|IU|) on very small scales indicates 
that a Lyapunov exponent will not be well defined. Clearly, this question deserves further 
investigation and we do believe that it is of considerable interest to understand in which 
sense we can or cannot have a "continuum limit" of a spectrum of Lyapunov exponents for 
a (non-dissipative) field theory like the SU{2) Yang-Mills theory when we study, say, the 
dynamical behavior in a given finite (physical) volume ~ L'^ ~ (Na)^ with periodic bound- 
ary conditions in the limit of long time t oo. A question of equal interest which may 
be related to the question above is to understand the (physical) spatial correlations in the 
system. Studies of (semi) classical dynamically generated temporal and spatial chaos may be 
important for the understanding of randomness in space and time of field configurations in 



Q.C.D. (see also discussions in e.g. [T^, |2g]). The present investigation has been devoted to a 



discussion of a possible "continuum limit" of (chaotic aspects of) purely c/asszca/ Yang- Mills 
fields regularized on a spatial lattice, whereas Yang-Mills fields in Nature are implemented 
as quantum fields, regularized by h (supplemented by a regularization/renormalization pre- 
scription). A connection between (semi) classical aspects of the real quantum theory and the 
behavior of the fields in the lattice regularized classical field theory would be most interesting 
to establish on a more rigorous basis. 
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Appendix 

We shall briefiy sketch the proof that for an SU{2) lattice Hamiltonian model (^ with a 
fixed lattice constant [a = 1), the maximal Lyapunov exponent is uniformly bounded from 
above, independent of the energy and the lattice size. We refer to for details in the proof. 

First we shall illustrate the main trick which is quite trivial: Consider a harmonic os- 
cillator whose evolution is governed by the equations x = y and y = —uj'^x. This system 
is integrable and thus it has vanishing Lyapunov exponents. Now, a way to show this is 
to consider the norm on the tangent space (which here is the same as the space itself) 
||(x, y)||^ = tu^x^ + y'^ which for uj non-zero is equivalent to any other norm on R^. The 
time-derivative of this norm is given by ^IKa;, ?/)|p = luj^xy + 2(— cu^x)?/ = and hence the 
norm is invariant under the fiow. It follows that all Lyapunov exponents vanish. 

In the case of SU{2) lattice gauge dynamics, the momentum P will play the role of u 
above and thus be part of an integrable dynamics (physically one may think of a rotor on 
a sphere). From this integrable part, we shall construct a 'good' metric for which norms of 
vectors in the tangent space would be constant in time, had it not been for the non-linear 
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coupling through the plaquettes. Taking the time derivative of the norm, we get a variety 
of terms which, due to the compactness of SU{2) and hence uniform boundedness of the 
couphngs, can be bounded by our 'good' metric itself ! This kind of bootstrap argument will 
give us a uniform upper bound as to how fast the norm can grow, and from this we deduce 
rigorous bounds on the maximal Lyapunov exponent. 

The Lie group SU{2) will be considered in the quaternion representation, which is a 
four dimensional embedding using Pauli matrices (not explicitly shown here) for which one 
writes : 

and define the usual Euclidean inner product in by : 

{U, V) = UqVq + UiVi + U2V2 + U3V3 = ^tT{UV^) . 

This yields a representation of the 6 dimensional tangent manifold T SU (2), such that a base 
point U G SU{2) and a tangent vector P G TuSU{2) are given coordinates (f/, P) G R'^Q)R'^, 
subjected to the constraints : 

{U, U) = 1 and (P, U)=0 . 
The entire phase space is denoted 

M = Y[T SU{2) . 

A 

Using the above notation, the Hamiltonian can be written as follows : 

i/ = Ei(P„P,)+^(l-itrf/o) , 

where the lattice constant is set to a = 1. The variation of H with respect to Ui gives us 
SH = {SUi, Vi) and thus a force acting upon Ui : 

The sum above extends over all neighboring links, such that (ijkl) forms a plaquette. Hence 
if the dimension of the space is d, Vi depends in total of 2((i — 1)(4 — 1) = 6{d — l) neighboring 
links but not on Ui itself. 

Under the constraints on U and P given above, the Hamiltonian gives rise to a flow 
(f)^ : M ^ M, t & R determined by the following differential equations (see e.g. : 

if- = p 

P, = V,~U,{U,Vi)-U{Pi,P,) . ^ ' 
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Consider a small deviation (5X of a trajectory 

X + 5X = {Ui + 6Ui, Pi + 5P,}ieA e M 

and linearize the above equation (0) in 5X to obtain the evolution now for a tangent vector 
which we again, by slight abuse of notation, denote (5X = {xi,yi}i^A ^ T^M 

yi = Vi-Ui{Ui,Vi) - Xi{Ui,Vi) 

- f/,(y„ Xi) - Xi{Pi, P,) - 2Ui{Pi, Vi) . 

A 'good' norm is now given by the following Riemann metric : 

^(^X, 51C) = W X,) (C + (P„ Pi)) + (y,, y,) - {x„ Pif - {y„ U^f) . (13) 

with an optimal choice of C = 2y/T7(d — 1). For any fixed energy, H, {Pi, Pi) is bounded 
from above by 2H, and a small calculation shows that G is positive definite as a quadratic 
form on TM. As we have noted earlier, our phase space is compact for any fixed energy and 
it follows that \fG provides a norm which is equivalent to any other norm on TM. If ||5X|| 
is another norm, compactness implies the existence of A;2, A;i > such that : 



A;i||5X|| < ^G'((5X,5X) < A;2||5X|| . (14) 

In particular, the maximal Lyapunov exponent measured with any two equivalent norms are 
identical. From ( |T^ ) we have : 

lin, 1 log - hm 1 log ^((^^'^^)[^]) (15) 

/™ t ||5X[0]|| " 2t G((5X,5X)[0]) ■ ^^^^ 

Now, using the constraints above, one verifies that the derivative of G along the flow equals : 

^G((5X,5X) = ^(x„y,)(C-(f/„V,)) + (x„x,)(P„V,) 

- (xj. Pi) {Vi, Xi) + {vi, yi) - {vi, Ui) {Ui, yi) . 

The Vi entering the equations above is the tangent vector associated with Vi and hence, it is a 
sum of 6{d — 1) terms of the type XjUlu} involving the tangent vectors Xj of the neighboring 
lattice sites (cf. above). A tedious but straight-forward calculation shows that the above 
time derivative satisfies : 

|^G(5X,(5X)| < 2AijG((5X, 5X) 



with Xb = y/d--l\/l + \/l + j^. One thus has the bound on the growth rate of the norm : 



G((5X[t],(5X[t]) < exp(ABt)^G(5X[0],5X[0]) . 
13 



In the case of c? = 3, we have = 4.03... which therefore gives an upper bound to the 
maximal Lyapunov exponent, independent of the energy and the lattice size. 

The metric we have constructed above is not gauge invariant and in principle one would 
like to study the dynamics and the metric of the associated moduli space, i.e. the space 
M modulo gauge transformations. However, in the study of upper bounds for Lyapunov 
exponents, this distinction is not necessary since gauge transformations are symmetry trans- 
formations of the dynamics and hence the associated Lyapunov exponents vanish. This is 



nicely illustrated by Gong in |T6[, where the complete Lyapunov spectrum is calculated for 
2^ and 3^ lattices. One third of the Lyapunov exponents vanish due to the gauge symmetry, 
expressed by the Gauss law constraint and residual static gauge transformations. A com- 
pactness argument again shows that using a gauge invariant distance measure (like the one 



used in |T1, ITS) which is only semi-positive but which distinguishes gauge non-equivalent 



configurations will yield the same measurements of Lyapunov exponents as a strictly positive 
metric. 
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Figure captions 



Figure 1(a) shows for SU{2) a plot of the maximal Lyapunov exponent as a function of the 
average energy per plaquette. The data points (diamonds) are adapted from Miiller et al. 
[p^ (see also Biro et al. [0, p. 192) where results are obtained from a simulation on a 20^ 
lattice with periodic boundary conditions in all spatial directions. The solid line is a linear 
fit through the origin and the dotted line is the function | x 0.38E^^'^ which is half of the 
result (cf. formula ([TTl) ) obtained from the homogeneous case (see also figure 1 (b)). 

Figure 1(b) displays a log-log plot of the same data points as in figure 1(a). On the figure 
we also show the (theoretically obtained) upper bound for the maximal Lyapunov exponent. 

Figure 2(a) shows a plot similar to figure 1(a) for U{1). The data points are adapted from 
Biro et al. [0. We note, that the Lyapunov exponent grows, as a lattice artifact, in exactly 
the same region of energy/plaquette 1 < E < 4 where the figure 1(a) (or 1(b)) is reported 
to exhibit continuum physics. 
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